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ABSTRACT 


The background of microstructure and its effect on sound propagation 
is reviewed. The expendable bathythermograph (XBT) system and current 
digitizing procedures are evaluated to determine their adequacy in 
investigating vertical thermal microstructure. Several methods of 
microthermal analysis are proposed. Temperature gradient versus depth 
plots are used to analyze data from several ocean regions - Atlantic 
Coastal, Andaman Sea, Eastern North Pacific, and the northern boundary 
of the Gulf Stream. Analysis of these data suggest a water mass depend- 
ence of the thermal microstructure. The significance of the small-scale 


structure to Naval operations is discussed. 
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I. SIGNIFICANCE OF MICROTHERMAL STRUCTURE 


Naval interest in the thermal structure of the ocean is presently 
centered on gross features such as sea surface temperature, mixed layer 
depth, thermocline gradient, and depth to the bottom of the thermocline. 
Some attention has been given, especially in the Canadian Navy, to the 
transients in the mixed layer [Tully, 1964] and their possible associ- 
ation with internal waves. The above gross features are used primarily 
as a means of classifying thermal structure. 

Technological developments have led to measurements which indicate 
the thermal structure is made-up of layers of varying thickness; from 
less than a centimeter to tens of meters [Stommel and Federov, 1967; 
Osborn, 1969; for example]. 

Of the variables affecting sound velocity in the ocean, temperature 
has the greatest effect. The presence of variations in temperature, or 
temperature gradients, causes sound velocity fluctuations to occur. 
Research on the effects of the thermal layers on sound in the ocean show 


that further investigation in theory and experiment is required. Small 


scale sound velocity variations cause irregular bending of sound rays. 





_ diverge the energy in a sound beam [Meyer and Romberg, 1963]. The 


physical features of the microthermal structure must be adequately 
described before their effects upon acoustic propagation can be fully 


understood. 


The purpose of this thesis is to review the present knowledge con- 
cerning microstructure as measured by expendable bathythermographs (XBT), 
to evaluate the capability of the XBT in measuring microthermal structure, 


to develop a scheme of quantitatively examining the microthermal structure, 


and to test this scheme on real data. 


II. REVIEW OF PRESENT KNOWLEDGE 


A. MICROTHERMAL STRUCTURE IN THE OCEAN 

Prior to 1937, vertical thermal structure was measured at discrete 
points using reversing thermometers. This method provided a point-to- 
point analysis of the vertical structure at one location. Even using 
point-to-point measurements, it was clear that the temperature vari- 
ability decreased with depth in most regions. In fact, the most common 
analysis of the temperature and salinity data, T-S correlation, usually 
ignored the upper 100 meters or so. However, this same analysis showed 
that in many ocean regions, in deep water, (e.g., below 1000 to 2000 m), 
the variability was on the order of the instrument error [Sverdrup, 
Johnson and Flemming, 1942]. This result was not unexpected because 
the gradients in the deep ocean became smaller. The discrete measure- 
ments of the Nansen cast are not suitable for studying small scale 
variability. 

The development of the mechanical bathythermograph (MBT) by Spilhaus 
[1938] provided the first convenient means of measuring a continuous 
profile of temperature. The mechanical bathythermograph was widely used 
during World War II. A large number of bathythermograms with double 
traces, which could not be explained as instrument errors, were dis- 
covered [Urick and Searfoss, 1948]. The double traces indicated that 
small-scale horizontal and/or vertical temperature variations occurred 
between the lowering and raising of the MBT. Urick and Searfoss [1949] 
recognized that this small-scale microstructure could be significant to 


sound propagation, depending on the size and shape of the inhomogeneity. 


Horizontal microthermal structure was measured by Holder in 1944 
using a thermopile mounted on a submarine. He detected horizontal 
temperature differences of 0.02F over distances of about 10 m [Smith, 
1967]. Measurements of these horizontal temperature variations were 
also made by Urick and Searfoss [1948, 1949] in July of 1948 near Key 
West, Florida. With a thermocouple mounted on the conning tower and, 
later, its bow, a submarine proceeded at a constant speed and depth, 
while the probe measured the water temperature to a resolution of 


0.Q03C. 
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Figure 1. Microthermal Variations at Various Depths 
[Urick and Searfoss, 1948]. 


Urick and Searfoss [1949] calculated the peak-to-peak root-mean- 


square (rms) variation for peaks in excess of 20 percent of the total 


10 


variation in a record calling this the root-mean-square maximum (rmsm). 
A significant difference between the value of rmsm in the mixed layer 
and in the thermocline was observed. In the mixed layer a mean rmsm 
value of 0.013C was observed. In the thermocline, the mean value was 
0.450C. The rmsm increased with depth in the mixed layer and decreased 
with depth in the thermocline. Applying Taylor's hypothesis, the mean 
horizontal thermal "patch size'' was found to vary from 5 to 30 m in the 
mixed layer and from 30 to 200 m in the thermocline region [Urick and 
Searfoss, 1949]. 

Liebermann [1951] mounted a sensitive platinum-resistance thermom- 
eter and a thermocouple on the periscope of a submarine. The platinun- 
resistance element could resolve spatial temperature variations on the 
order of 10 cm in horizontal distances as the submarine proceeded at a 
constant speed of a few meters per second at a fixed depth (Figure 2). 
To characterize the temperature fluctuations, he used the autocorrelation 
function. "Patch size" was defined as the distance of the microthermal 
structure, within which the temperature variations maintained coherence. 


The function is given by: 


oe 


R (0) = [ [7 (x): Getog)] ax. (1) 
0 


His results could be fit equally well by an exponential, 


-op/a 
ReG@e="e a (2) 
or the Gaussian, 


2 
aig eecmeae (3) 


where: 


Pt 


T (x) temperatute at some point x; 


z (xto 9) = temperature at another point displaced a distance Po 
Peon xX; 
© = autocorrelation distance; 
<p> =a = "patch size". 


Whereas, the exponential function fitted the data points and decayed 
rapidly as p increased (Figure 2), it did not have a zero slope at the 
origin, thus, the exponential function was not appropriate for very 
small correlation distance (e.g., a few cm or less) because of heat 
conductivity in sea water. He proposed the Gaussian function to 
correct this difficulty, reasoning that thermal conductivity would 
prevent any discontinuity in the temperature field [Liebermann, 1951]. 
Liebermann assumed that the thermal microstructure represented a 
cloud of spheres with a Rayleigh distribution having a mean size value 


of the correlation distance and therefore yielding 


R (0) =R (a) =e. (4) 


Using this definition, the mean size of the inhomogeneities, a, was 
found to equal 60 cm in the ocean (Figure 2). 

Piip [1961] investigated the sound velocity microstructure in the 
sound channel off Bermuda with a NBS velocimeter. It was noted that a 
number of thin steps remained over periods of hours, some oscillating 
30 m in depth in a three-hour period. At the same station and same 
depth, sound velocity changes of one m/sec in a few hours were common 
and "nearly the rule'’. The sound velocity axis was observed to change 
150 m in depth and 3.5 m/sec over a three-month period. Pipp showed 


that the base of the main thermocline (about 1000 to 1300 m) exhibited 


eZ 


relatively large and rapid changes in sound velocity in an inhomogeneous 


region in a constant state of fluctuation. 
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Figure 2. Temperature Fluctuations and Associated Correlation Function 


(After Liebermann) 


Small-scale temperature and salinity structure was measured with a 
salinity, temperature, and depth recorder (STD) by Stommel and Federov 
[1967] near Timor and Mindinao during July of 1965. Their measurements 


revealed that the vertical microstructure actually consisted of very 


The thin layers were separated by thicker 





layers of nearly homogeneous (isothermal or isohaline) water, i.e., 


step-like in appearance on their STD traces. By following a discrete 


13 


feature from one cast to the nexti|homogenéous layers were found to be 





‘two to 40 m (Figure 3). A boot-like feature was found at the bottom of 


the mixed layer (Figure 4). At a depth of 200 to 300 m, the mean thermal 
gradient was about 0.040C/m. In the thin steps, the gradient was less 
than 0.003C/m, and in the homogeneous steps, the thermal gradient was 
greater than 0.350C/m. Stommel and Federov concluded that the ocean 
consisted of nearly homogeneous layers separated by steps with extremely 
high thermal gradients. 

Tate and Howe [1968], in the summer of 1966 in the Northwestern 
Atlantic, used an STD to look at the thermocline stratification in the 
Mediterranean Water intrusion at a depth of 1280 to 1500 m. They found 


ture, from 0.02 to 0.055 o/oo in salinity, and from thickness of meme 


30 m (Figure 5). SAnatep's°@¥@tag@Nt maa Miiiaashangesv RUSS 


ao, change of -0.044, while the average salinity step of 0.044 o/oo 


produced a 0, change of +0.039. The resulting change of O, (-0.005) 
gave a slightly stable layer system in a situation which was clearly 
one of opposing salinity and temperature density gradients. 

Cooper and Stommel [1968], using an STD, found the main thermocline 
in the Sargasso Sea to consist of regular salinity and temperature steps 
of homogeneous layers three to five meters thick, separated by trans- 
itional layers 10 to 15 meters thick, with temperature changes of 0.30 
to 0.50C and salinity changes of 0.04 to 0.10 o/oo (Figure 6). The 


horizontal extent of these steps was between 400 to 1000 meters. The 
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Figure 3.a. 


Figure 3.b. 
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a. b. 


Photographic copy of actual STD trace (not retouched) of 
station 21 off Mindinao. Depth scale is 500 m, each small 
division being 5 m. The salinity curve, labeled S, has 
the smallest division being 0.05 o/oo. The temperature, 
labeled T, has a smallest division of 0.05C. The lamina 
B is so labeled. Because the pens cannot both traverse 
the same path mechanically, the temperature curve is 
offset downward by 5 m and always reads 5 m too deep. 
The salinity scale is not offset in depth. Uniform 
corrections for salinity of -0.9 o/oo and temperature of 
0.07C have not been applied. 


A composite chart-graph for lamina B off Mindinao. Each 
square is 0.05C wide (8.5 to 9.0C) and 50 m high. The 
small arrow indicates a depth of 400 m in each square. 
Numbers in the lower left hand corners are the station 
numbers. The distances from stations 18 to 25 and 26 to 


35 are about 10 nm. 


[Stommel and Federov, 1967] 
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Figure 4. From the soundings in the mixed layer off Timor, gross dif- 


ferences in the properties in different locations are self- 
evident. There is considerable variation in vertical 
structure but each sounding shows a nearly homogeneous 
region near the surface (station 13) but some stations 


(like 7) have shallow irregular layers overlying the mixed 
layer. 


[Stommel and Federov, 1967 ] 
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Figure 5. A typical step in Mediterranean Water [Tait and Howe, 1968]. 
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Woods [1968, 1968a] investigated the microstructure in the summer 


thermocline off Malta using a temperature-gradient meter for obtaining 


the thermal gradient between tw 





Gilalvertical staff. The resolution of the thermal gradient was within 


0O.01C per thermistor separation distance. An additional thermistor 
recorded temperature. Woods plotted profiles of the thermal gradient 
and temperature versus depth (Figures 7, 8 and 9). Using time-lapse 
photography, Woods obtained pictures of streaks left in the wakes of 
free-dropping dye pellets for shear measurements. A dye-packet array 
tied between moored, submerged floats injecting dye into selected levels 


of the thermocline, indicated that the ae 





only centimeters thick. (Sometimes, he detected sheets up to a meter 


thick but they appeared as aggragates of the 10 cm thick sheets.) The 


WRN oc oscanbalence.. Mature sheets contained thin laminar-flow 


zones in their centers. Neighboring layers had dissimilar properties 
indicating different salinities, velocities, and turbidities. This 


suggested different trajectories and a very slow exchange of properties 
across the intervening sheets. If turbulence was small, 


flux exchange 
was small; therefore, the exchange had to rely on molecular processes, 


which were very slow indeed. 


Individual layers were followed horizontally for distances up to 


72 km by soundings from high-speed launches. Layers were identified on 


Ey 


(uz) yydeq 


"(8961 ‘Teumo0js 


‘[8961 ‘Spoom] Sutoeds 10IsTurayQ pue 1sedo009] days uSNOsuegouoy,, ay Burs 
m2 OS & YATM spoom fq peutejgo B SsO1de UOTSsad.qns ptde1z ut s3utpunos 
SON] eA yjdap-quat pea3 291Injerisduaj Ajtutyes pue evinjeisduiaq [erzaaas SutyeQ 
pue yjdop-sanjeaaduay Teotddy y -y7 aan8ty Aq apeu eas OsseZeieg ay} ut sdaqc "9 21nB1Ty 








WOZ9 Yydap aboisay 


(aro 9¢/5,) quoIpe12 einqyeieduiay, 








ie 
he ee 
rage 
i == 
tee 
fe 
o 
ee 
| a xr. 


0£°90 99 "MHA “4z x4 Ss Deal 


18 


Temperature gradient 










(deg C/25em) Temperature (°C) 
02 04 66 20 2! 92 93 4 25 2 2) 
Wave-mixed 
layer 
re Tabet 
Surface 
layer 


Se ewseseeees eee weeenecun Seoacesavas 


20 
—F eee 
3 v7, 
= 
© 
& f 
€ 
& 
30 
Se hermocline 
/ lay ers 
50 
Figure 8. for Maltese 


Waters, 14 September 1967 [Woods, 1968a]. 
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Dye was spread 
in this region by internal waves to form a "broad carpet of color". This 
colored sheet made it possible to photograph the internal waves as they 
crossed it. Wave-lengths between 5 and 10 m were dominant and appeared 
in long, coherent trains, each associated with a single sheet. They had 
heights up to 1.0 m, a phase speed of 2.5 cm/sec and periods of about 

5 min. These waves only affected a single sheet, being much shorter than 
the long internal waves which move the entire thermocline. 

A sheet was defined as stable when it appeared smooth everywhere ; 
this indicated laminar flow. When a wave passed over a sheet, shear had 
its greatest effects at the crest and trough. Wave shear and drift 
shear interacted either cancelling or reinforcing each other causing a 
sheet to become unstable and patchy. 

Wavelets, about 75 cm in length, caused banding to appear parallel 
to the crests and troughs. These wavelets generally formed, grew, and 
broke in a period of less than two minutes. Upon breaking in classical 
fashion, a second smaller breaker was generally thrown forward. The 


wavelets had a maximum height of about 20 cm. After breaking, a turbu- 


lence region remained for about 5 minutes, and the scar in the sheet 


remained for several hours. These ‘scars were a common feature in 





Za 


examples of sheets. Woods divides the thermocline into a half-dozen or 


so layers of low shear and moderate thermal gradient separated by thin 
liminar flow sheets on intense shear and thermal gradient. He discussed 
the important aspect of heat transfer through such a layered ocean. He 


suggested that the molecular heat transfer through the high gradient 


stable layers was negligible, but that heat transfer does take place 


The use of the Ramsey SVTP (sound velocity, temperature, and pressure) 
probe provides a digital output of these parameters. This instrument has 
an accuracy of + 0.01C and + 0.25 percent of depth with a precision of 
+ 0.08 m/sec in sound velocity [Lovett, 1968]. In view of the plus-or- 
minus-one count for electronic counters, small-scale features, i.e., 
less than + 0.01C and + 1m in depth, are ambiguous. Lovett [1968] 
discussed the thermal gradients determined with this instrument (Figure 
10). The 10 m depth increments he used eliminated some microstructure; 
however, using 5 m increments introduced random readout errors. 

Ingham [1968] used an STD probe to investigate the mixed layer in 
tropical waters. The vertical salinity and temperature structures 
forming the "isothermal" and "isohaline'' layers were generally analogous 
in the near-surface layers. The difference in the depths of these two 
layers exceeded 7.5 m in 45 percent of the casts and 17.5 m in 31 percent 
as shown by Figure 11 [Ingham, 1968]. Microstructure effects on sound 
velocity in the near-surface layers would depend on salinity in some ocean 


regions. 
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Figure 10. Temperature Gradient per 10 m, Mean Gradient per 100 m, and 
Average Deviation from Mean versus Depth from Ramsey Probe 
Data [Lovett, 1968]. 
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Figure 11. Histogram of Depth Differences AZ of the Isothermal and 
Isohaline Layer Depths [Ingham, 1968]. 


23 


Pingree [1969] using an STD, recorded analogue signals on tape. He 
made fine-scale plots of salinity and temperature. In particular, he 
investigated the stability of the deep water, i.e., about 500 to 2000 m 
(Figure 12). 

“Osborn [1969] using a free-falling temperature gradient meter 
ta ce Ete ent (J Off 
San Diego, the horizontal extent of the layers was greater than 750 
meters in the seasonal thermocline, but only a few hundred meters at 
depths greater than 400 meters. He noted the layers were, in fact, 
made up of even finer layers with a thickness of only a few centimeters. 
His average sheet was about one meter thick. 

Recent work by Neshbya, Neal, and Denner [1969] using the expendable 
bathythermograph manufactured by General Motors Defense Laboratories, 
obtained records of temperature steps (Figures 14, 15 and 16). While 
on Ice Station T-3 in the Arctic (near 84-38N, 128-22W), they lowered 
XBT probes by hand. The Ice Station results are the first XBT micro- 


structure records without ship motion. 


B. APPLICATIONS 

Urick and Searfoss [1949] calculated theoretically that a sound ray 
would be bent through an angle of only 3.7 sec by the average inhomo- 
geneity in the mixed layer but as much as 2.1 min of arc in the 
thermocline. 

Sheehy [1950] measured sonar signal fluctuation as functions of depth 
and range. He found the coefficient of variation, defined as the stand- 
ard deviation of the amplitudes expressed as a percentage of the mean 


amplitude, varied as the square root of range and did not appear to be a 


function of depth. 
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Figure 12. T-S Curve from 1500 to 1800 m with a Value Every Meter. 


Circles Represent Increments of 2 m in Depth [Pingree, 
1969]. ; 
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Figure 14. Vertical Profile of Temperature Under Ice Station T-3 


(84-38N, 128-22W) 19 March 1969. [Neshbya, Neal and 
Denner, 1969]. 
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Time Series from 21452 17 March to 20332 18 March 1969. 
(Neshbya 1969]. 


Figure 15. 


Neal and Denner, 
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Figure 16. Time Series with High Amplification Taken from T-3 in March 
1969. Temperature increases from right to left. [Neshbya, 
Neal and Denner, 1969]. 
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“acoustic propagation He wrote the sharpness of temperature boundaries 


determines influence of the microstructure on acoustic propagation. 


“energy depending upon the size and shape of the "patch'' The tempera- 
ture microstructure may reflect, refract, scatter, or focus acoustic 
waves, The microthermal structure was considered to be made up of 
randomly sized ''spheres'', spaced randomly throughout the medium, with 

a scattering cross section dependent upon the autocorrelation of the 
temperature inhomogeneities, R (9). Comparing derived scattering cross 
section with actual scattering in the ocean (i.e., volume reverberation), 
Liebermann concluded the observed reverberation was higher than one 

would expect from thermal inhomogeneities He attributed volume 
reverberation to biological scatters. 

Liebermann found that refraction effects of the inhomogeneities 
resulted in warping of thé wave fronts, and caused the wave front to 
be a complex surface exhibiting both convergence and divergence. The 
temporal motion of the inhomogeneities resulted in sound "scintillation" 
at the receiver. 

Mintzer [1953] found Sheehy's coefficient of variation to be depend- 
ent upon the three-halves power of the range rather than the one-half 
power of range. 

Skudrzyk [1957] discussed the focusing and defocusing effect within 
the focusing range Ty given by r, = tae where a is the "patch" radius 


0 


and k is the acoustic wavenumber At ranges greater than Ty» in the 


interference range, phase interference results (Figures 17 and .18). 
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Figure 17. Illustration of the Various Types of Acoustic Scattering 
Due to a Thermal Patch in the Sea [Whitemarsh, Skudrzyk 
and Urick, 1957]. 
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Figure 18. The Phase Angle Between the Incident and the Scattered 
Pressure and the Physical Nature of the Scattering as a 
Function of the Range [Skudrzyk, 1957]. 
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Using this concept, Skudrzyk (Figure 19) shows the focusing range Tr) as 
a function of frequency and patch size. (For example, a 50 kHz frequency 


Sonar has a focusing range of 75 meters, if the patch size is assumed to 


be 60 centimeters.) For finding patch a 


Further, Skudrzyk started theoretical investigation with spheroidal 





patches, which he found to have a considerably stronger focusing effect 
than the spherical patches. 

Whitemarsh, Skudrzyk, and Urick [1957] did additional work with the 
Spheriodal theory. The correlation functions previously used were dis- 
cussed. It was concluded that the present correlation functions were 


not suitable and either had to be replaced or redefined. The authors 


argued if the temperature microstructure was produced by turbulent 
This, they show, was in excellent agreement with their data. 


lending support to continued studies such as are presented in this 


thesis, their final conclusion states that the 
a greatest gap seems to lie in the microthermal 
Statistics of the ocean. If these were known better 
than they are now, apparently the amplitude statistics 
of the fluctuations of sound transmitted through the 1 
sea could be predicted roughly, but with some assurance." 


Works by Peterson [1965], Barakus [1968] and DiNapoli [1969] have 
considered the normal mode approach to the interference range solution 
for the temperature inhomogeneities and have also investigated the 


effects of internal waves on sound in stratified mediums. 


Nan venueaane D. C., Skudrzyk, E., and Urick, R. J., "Forward Scattering 
of Sound in the Sea and Its Correlation 7 the Teieer fears Micro- 
Structure,'' J. Acous. Soc. Am., v. 29, p. 1141, October 1957. 
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C. INADEQUACIES OF THE PRESENT THEORY 

Since the thermal microstructure is not related simply to the gross 
thermal features measured on MBTs, there is presently no easy way to 
predict, on a day-to-day basis, the effects of microthermal structure 
on the propagation of sound in the sea. The theories, so far, have 
developed in two classes. 

The first theories [Liebermann, 1951; Mintzer, 1953, 1953a, 1954; 
and Skudrzyk, 1957] concerned spherical and, more recently, spheriodal 
shaped inhomogeneities throughout the ocean. Part A above has dis- 
cussed the present knowledge of microstructure showing that the long 
thin sheets and layers do exist as microthermal structure. These first 
acoustical theories do not apply to what is now described as the micro- 
structure of the ocean. 

The theoretical treatment by Whitemarsh, Skudrzyk and Urick [1957], 
considering speriodal patches and Kolomogorov's turbulence theories, 
could not find any real theoretical explanation of what happens in the 
sea when sound passes through an inhomogeneity. The investigations 
using normal modes have not produced any practical operational results. 

lee [1961] using ray path theory, demonstrated the effect of a 
typical sinusiodal internal wave on sound intensity in a three-layered 
model. For the medium with an internal wave, intensity variations up 
to 22 db were found in distances of less than one internal wavelength, 
as shown by Figures 21 and 22. In contrast, the same medium, without 
an internal wave, produced a 5 db variation, including a 2.5 db 
spreading loss, over the same distance. 

Another study investigation the effects of sound velocity micro- 


structure was conducted by FNWC [Wolff, Tatro and Megehee, 1967]. A 
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Figure 21. Diagram of Sound Rays through a Medium Which has an 


Internal Wave on the Thermocline [Lee, 1961}. 
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Figure 22. Sound Level in the Medium with an Internal Wave on the 


Thermocline. The db reference level is that corresponding 
to a sound level of 60 db at 1 foot from the directional 
source (X equals 0, Z equals 10 ft) along the horizontal. 
The field is contoured at intervale of 2.5 db [lee, 1961]. 
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sophisticated ray-tracing program [Ayers, Wolff, Carstensen and Ayers, 1966] 
was used for the analysis. Three complete ray trace runs were made to 
the first convergence zone using 120, 54, and 20 digitized sound velo- 
city points from one SVTP trace containing sound velocity microstructure. 
The 120-point sound velocity profile preserved all the "wiggles" or 
microstructure, the 54-point profile used values at every 50 meter of 
depth, and the 20-point profile used the standard oceanographic depth 
values for sound velocity (Figure 23). The range to the first conver- 
gence zone agreed within 150 yards for all three cases. Based on the 

ray trace runs only, it was concluded the sound velocity microstructure 
was "essentially self-cancelling'" and did not affect the convergence 

zone parameters {Wolff, Tatro, and Meghee, 1967]. Propagation loss and 


reverberation effects were not reported. 
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Figure 23. A Portion of the Three Different Sound Velocity Profiles 
Used by FNWC for Ray Tracing [after Wolff, et al.]}. 
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III. MEASURING MICROTHERMAL STRUCTURE 


A, THE BATHYTHERMOGRA PH 
1, Mechanical Bathythermograph 

At the present time temperature is the most readily observed 
variable in the upper few hundred meters of the ocean. A major con- 
tributor to this situation is a singularly successful oceanographic 
instrument known as the mechanical bathythermograph or MBT. Subsequent 
to its introduction by Rossby and Montgomery [1934] and perfection by 
Spilhaus [1938], the MBT has become one of the basic instruments of 
almost all oceanographic agencies because of its simplicity of 
Operation, reliability and relatively low cost. The result has been 
that a large number of bathythermograms have been collected. 

Researchers during World War II, using the MBT, observed micro- 
thermal structure on their traces. No theory had been formulated which 
took into account the effects caused by microthermal structure on under- 
water sound transmission. It was recognized, however, that temperature 
microstructure must be considered in any explanation of observed under- 
water sound transmission [National Defense Research Committee, 1946]. 

2. The Expendable Bathythermograph 
a. Capabilities 
Presently the U. S. Navy is using the XBT system manufactured 
by the Sippican Corporation of Marion, Massachusetts. This system con- 
sists of an expendable shipboard cannister and probe, shipboard launcher, 
and a strip chart recorder. Data so obtained are submitted to FNWC and 


were used in this thesis. 


oe 


The original Bureau of Ships' specifications [{Denner, 1966] 


for an XBT system are given in Table I. 


TABLE [I 

Temperature range 28F to 96F 

Temperature accuracy + 0.4F 

Depth range O to L0C0 fe 

Depth accuracy hi 2% or 15 ft, whichever 
is greater 

Ship's speed when launching ae 63% of a temperature 
step after 3 ft depth 
change 

The expendable Sippican probe currently used has a 
thermistor sensor in its nose. The probe is connected electrically 
to the recorder by two conductor wires which pay out from both the 
probe and the canister. The cannister in the launcher is connected 
to the recorder by a cable. There are no splices in the wire. The 
probe's weight is a critical factor in determining its fall-rate. 

The Sippican system was found to be within + 3 gms both with and 
without its internal copper wire. (The deviation was from 0.3% to 
0.64% of the total weight of the probe.) The observed average fall- 
velocity of the probe was 20.3 + 0.2 ft/sec. A probe, with its 
internal wire, had an average terminal velocity of 19.9 + 0.7 ft/sec. 
Terminal velocity was reached at 18.0 + 0.5 ft [Gouzie, Sanders and 
Littlehale, 1966]. 

The present XBT strip chart recorder is a Mark 2A model 
using pressure sensitive paper for recording the thermal profile. The 
paper moves vertically on the recorder at a rate of 3.25 in/min. In one 
minute, the XBT, falling at 20 ft/sec, will have reached a depth of 1200 


ft. Thus, 1200 ft will be recorded on 3.25 inches of chart papergeon 


0.0027 of an inch of chart paper is used for every foot of fall. 
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The chart paper records horizontally from 28F to 96F in 
7.0 inches, or about 0.103 of an inch of chart paper per 1.0F. Table 
II shows the resolution capabilities of the paper using different 


recording intervals. 


TABLE II 
Resolution 
Chart Spacin inches Temperature (F Depth (ft) 
0.01 0.0971 3.690 
0.005 0.0476 i045 
0.001 0.0097 .369 


The response time for the thermistor is 110 milliseconds 
nominal and 130 milliseconds maximum [Choate, 1970 and Demeo, 1968]. 
A time constant-temperature change plot shows that 98.2 percent of a 
temperature change has been recorded in 10.4 feet of fall, 83 percent 
recorded in 5.2 feet of fall, 68 percent in two feet of fall and 35 
percent in one foot of fall (Figure 24). 

‘The thermistor, because of its finite size, has thermal 
inertia which causes it to act as a low pass filter. As a result, 
fine structure is lost or at best smoothed out. The original trace 
could be reconstructed if the filtering characteristics (of both phase 
and amplitude) of the thermistor were known. These could be obtained 
experimentally. 

b. Data Handling 

Present procedures for handling XBT data are indicated by 
the log sheets which accompany the XBT‘'s in the shipping container. 
First, a BATHY or BAXBT message is required to be submitted to the 


appropriate addressee using a PRIORITY precedence message [National 
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Oceanographic Data Center, 1968, 1968a]. Upon receipt of the BATHY or 
BAXBT message the receiving agency uses this as synoptic data for sub- 
jective ocean forecasts in his locality. The BATHY or BAXBT is for- 
warded through the Naval Environmental Data Network to FNWC in Monterey. 
FNWC uses the message reports to update its ocean fields for objective 
forecasting in various oceanographic and meteorological products. 

Upon completion of a cruise the ship which took the XBT 
cast mails the original trace and log sheet to FNWC for digitizing. 
The original XBT trace is digitized for use by FNWC, and, when 
requested, the results of the digitizing are provided to the XBT 
submitting activity. Upon completion of digitizing, the original XBT 
traces are forwarded to the National Oceanographic Data Center (NODC) 
for storage in their archives. NODC also receives a copy of the digital 
data on tape. 

Present procedures at FNWC use DIGITBT and several other 
specialized XBT data handling programs. New XBT digitized data are 
stacked on the XBT Master File (presently at a rate of about 500 traces 


per week). 


B. OTHER SENSORS 

sensors for measuring microthermal structure, other than those 
discussed elsewhere in this thesis, are numerous. 

From a fixed platform, like the NEL Oceanographic Tower, thermistor 
strings and arrays are used. Bottom devices, such as thermistors, 
attached to the shore or a tower by a wire, have been employed to give 


continuous temperature records. 
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From submarines, thermopiles, thermocouples, thermistors, platinum 
resistance thermometers, and thin film thermometers have been employed. 
Almost any oceanographic electronic or mechanical thermometer could be 


used on a submarine. 


C. DIGITIZING PROCEDURES 

At the present time, all XBT traces are submitted to Fleet Numerical 
Weather Central (FNWC) in Monterey, California, the only digitizing 
center for these strip chart analogue traces [National Oceanographic 
Data Center, 1968, 1968a; and Bauer, 1969). 

At the Point Pinos Annex of FNWC, two Calma Model 480 flat-bed 
digital recorders are used to digitize XBT traces in 0.01 inch 
increments. These digital recorders are operated in a mode which makes 
the recorded intervals dependent on the tracing speed. The Calma 
output recorded on digital magnetic tape containing binary mode 
Characters at a density of 556 bits per inch (BPI) with even parity 
and variable length records. 

FNWC program DIGITBT, written by Lieutenant N. L. Perkins, USN, 
in FORTRAN IV for FNWC's CDC 6500 computer, is used currently for 
the reduction of data from XBT traces. FNWC tape input, tape output, 
and some other subroutines of DIGITBT are written in COMPASS language. 
In order to analyze microthermal structure, DIGITBT was extensively 
modified. The new program, referred to as MICROXBT, is included as 
Appendix A. The basic procedures for DIGITBT and MICROXBT programs 
are the same. 

DIGITBT and MICROXBT programs both convert the X and Y measure- 


ments of the digitizer into depth (X) versus temperature (Y) values. 
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Corrections are made for nonlinearities produced by the XBT recorder in 
both the depth and temperature directions. All shipboard XBT recorders 
are assumed to produce error-free traces. The 62.0F temperature line 
drawn on the trace when the probe is loaded is used as a baseline for 
digitizing. No correction is made if this baseline temperature line 
is not the same as the 62.0F line printed on the trace paper. An 
eye-interpolated value of the XBT-drawn baseline temperature on the 
printed temperature grid is noted as a part of the identification 
record. 

The rigid procedures for digitizing require that exactly one 
identification (ID) record and two digitized trace records for each 
XBT be recorded on the tape. The ID record information contains the 
data (day, month, year); time (hour and minute); latitude and longi- 
tude (degree, minute, hemisphere); baseline temperature (to nearest 
O.1C or O.1F depending on the grid type); printed-grid type; and maxi- 
mum depth for the probe (coded as 15 for 1500 ft XBT or 25 for the 
2500 ft XBT). Each trace record contains pairs of increments (AX, AY) 
read from the baseline to the true surface temperature, a point "flag" 
and then pairs of values (AX, AY) along the trace to a cut-off depth. 
Each trace is digitized twice. The complete procedure is described 
by Bauer [1969]. 

To process digitized XBT data using the DIGITBT program, the two 
trace records are compared point-by-point in a deviation test. If 
the X and Y paired values (i.e., depth and temperature) differ by more than 
0.035 inches, DIGITBT rejects the whole XBT record and the XBT trace must 


be completely re-digitized. If the XBT record passes the deviation test, 
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the two traces are averaged, point-by-point. The resultant trace is 
used for the rest of the program. 

In the following the term "inflection point" will be used to pre- 
serve consistency with the author of DIGITBT. Inflection points are 


‘criteria as significant points 


those points which meet the programs 
which are the programs output. 

Both DIGITBT and MICROXBT select the first and last observations 
as inflection points in the output list. The two programs differ at 
the inflection point choice stage. Other observations in DIGITBT are 
screened in groups, and, for each group, the single observation is 
found which is most different from a linear interpolation between the 
two end points of this group of observations. Each observation in a 
group is measured to find how out-of-line it is in temperature and in 
depth, and the smaller of these two is considered the out-of-line 
measure. All observations in the group are then compared to find the 
one with the largest out-of-line measure. If that value exceeds 
0.0150 for an inch, that observation is selected for output in DIGITBT; 
otherwise, the group expands to include the next value and the inflec- 
tion test is repeated until ansanileeelongpoljeeee source. 

The selection of inflection points was modified in MICROXBT to use 
a minimum of 0.0001 of an inch for the out-of-line measure. This pro- 
duces an output containing every inflection point, which is nonlinear 
between two consecutive points in the digitized record. Any observations 
which cannot be reconstructed by linear interpolation are not required 
and are removed from the output listing in both DIGITBT and MICROXBT. 

DIGITBT converts depth units from digitized inches to ft of depth 


by linear interpolation using a depth table and from ft to m by formula. 
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Temperature units are converted from digitized inches to F by a tempera- 
ture table and from F to C by formula. These depth and temperature 
tables were constructed to account for the nonlinearities between the 
recorder and the actual trace. Values of depth are rounded to the 
nearest m and temperature to the nearest 0.01C in DIGITBT. MICROXBT 
has eliminated the ft to m and F to C conversions. Depths are rounded 
to the nearest 0.1 ft and temperature to the nearest Q.OOIF. 

After the conversion from ft to the nearest m, DIGITBT removes, 
for a second time, any of the rounded-off values which are nonlinear 
between two neighboring points. If depth values appear as duplicates, 
one value is omitted in the output listing. In MICROXBT, there are no 
duplicate values because 0.1 ft in depth is well beyond the hundredth 
Siwan inch digitizing capability of the digitizer. Thus, by reading 
depths in tenths of a ft and temperature in thousandths of a degree F, 
MICROXBT looses no values from duplicate depths, and the possibility of 
introducing false linearities is reduced considerably. For these rea- 
sons, MICROXBT gives a much finer digitized record than DIGITBT. 

The DIGITBT output tape is written in binary for the CDC 6500; the 
MICROXBT output tape is written for the IBM 360/67 at the Naval Post- 
graduate School (NPS). Appendix C contains the job control cards 
(JCL) required for using a CDC 6500 output tape from MICROXBT with 
external BCD at a density of 556 BPI as in input to a separate micro- 


thermal analysis program at NPS. 
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IV. MODELS FOR EXAMINING MICROTHERMAL STRUCTURE 


All models and methods used in this thesis required the output tape 
from MICROXBT to be used as a direct input on the seven channel tape 
unit. The XBT ID record information was read first and if the recorded 
XBT was from a desired geographical region the data (depth and tempera- 
ture) were analyzed for microthermal structure at NPS. 

A "best" value for the temperature structure was found by using the 
depth inflection points from MICROXBT. For this "best" array, tempera- 


ture was linearily interpolated to a value for every foot of depth. The 


mixed layer depth (MLD) was determined, using inflection points from 
MICROXBT, as that depth immediately above the depth where the temperature 


“was 2.0F different from that at the surface [Naval Weather Service, 1967]. 





Several microstructure models were tried in order to perform an 


objective analysis of the oceanic microthermal structure. ‘The use of 





A. DERIVATIVE MODEL 
1. Mathematical Model 


The most significant feature of the thermal structure is the 


thermal gradient defined as the derivative of the temperature with 


respect to depth (dT/dZ). By plotting the gradient versus depth (Z), 


a parameter such as the thermal gradient is related to stability, 





Furthermore , 


entrophy generation, turbulent heat transport, and the Brunt-Vaisala 


frequency [Osborn, 1969]. 


2. Experimental Methods 


There are two methods for experimentally determining the thermal 
gradient presently in use. 
a. Electronic 
A resistance thermometer in an electronic bridge passing 
through water produces a signal which is proportional to temperature as 
a function of time. By differentiating this signal with respect to time 
(which is analogous to differentiating with respect to distance), the 
values of gradient versus distance the probe travels are found. 
This method was successfully employed by Grant, Moillet and 
Vogel [1968], in the horizontal direction for temperature and turbulence 
measurements. A thin-film resistance thermometer, mounted on the bow of 
a submarine, measured temperature versus horizontal distance traveled. 
The time derivative was taken to produce a direct plot of dT/dX, the 
horizontal thermal gradients. 
Osborn [1969] used a thermistor, capacitively coupled to 
two differentiating amplifiers, to measure the vertical temperature 
gradient with his free-falling gradiometers. His resulting output 
(Figure 13) shows a gradient versus depth record for the North Pacific 
off the San Diego Trough. 
b. Two-probe System 
“Temperature gradients may be found by using two vertically 
‘spaced thermistors. The temperature difference is a measure of the mean 
gradient between them. Such a gradient-meter was used by Woods [1968, 
1968a] who obtained gradient at various depths (Figures 7, 8, and 9). 
This method is limited to the resolving power of the thermistors and is 


a function of thermistor spacing. 
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3—  Numericat 


Bradient but uses the Finite differeneelappreximativon’ The derivative 


of the temperature at any depth T(Z.) with respect to depth Z is defined 





as; 


dz L268 bin ee 
i i 
which can be shown to equal: 


dt Syren 


dZ AZ—0 AZ 


a Lime) E(Z bane T(Z) 


AZ~0 AZ 
For this thesis, 4Z = 1 foot, which yields the finite difference approxi- 
mation for the vertical temperature derivative: 
S = or. (7) 
B. MISCELLANEOUS METHODS 
1. Statistical Bias and Variance 

For this method, a standard or smooth curve was required. 
Several standard curves were investigated. The variance and bias were 
calculated between the standard and the "best'’ curves. The standard was 
found by computing the temperature values for a series of specified depths. 

The first standard curve investigated used MLD and the FNWC 
significant levels (surface, 100, 200, 300, 400, 600, 800, and 1200 ft). 
The temperature values were found using the "best" temperature values at 
those depths. Next, an array of temperature values was found at foot 
increments by linear interpolation between the temperature values at the 


significant levels or the MLD. This was designated the FNWC significant 


level curve. Calculations of bias and variance did not reveal any cor- 
relation with the microthermal structure easily visible on the XBT trace. 

The second standard curve investigated used the temperatures at 
the standard oceanographic station depths (surface, 10, 20, 30, 50, 70, 
100, 150, 200, 250, 300, 400, 500, and 800 m). Since these depths are 
metric, a linear interpolation was done between the two neighboring 
values on the "best" curve. Values of temperature for every foot were 
found by linear interpolations between the standard depth temperature 
values. Again, no significant correlation was noted. 

The third standard curve used thse standard depth temperature 
values but used these standard depth values as an input to the LGINTP 
subroutine [Denner, 1969]. This subroutine uses the four standard depth 
points neighboring the depth in question to make an "average" curve from 
an upward and downward parabolic interpolation and a linear interpolation. 
This scheme produces a smoothing of points on the curve. This technique 
is usually used for interpolating standard depth oceanographic data 
[Rattray, 1965]. These values, after computing the bias and variance, 
again gave no correlation with the microthermal structure. 

The LGINTP subroutine was not used with the FNWC significant 
level values for smoothing since the FNWC fields assume linear inter- 
polation. 

Other standards considered, but not investigated, included a 
Gaussian standard, using the technique for approximating the temperature 
structure proposed by Grosfils [1968] and an exponential standard. 

2. Gradient Distribution Method 
For investigating the temperature structure, the "best" tempera- 


ture data were separated into the FNWC significant levels. For each 
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Significant level, the temperature change from one depth to the next 
were calculated and then sorted in numerical ascending order with the 
IBM subroutine SHSORT. With these values of temperature change per foot 
or gradient (F/ft), the Cemperature structure was now characterized by 
the individual thermal gradient values and a total temperature change AT. 
Because of the range of temperature values in the thermocline, this region 
of the XBT trace was chosen for developing this method. 

The capabilities of the XBT indicate a gradient interval of 0.02 
F/ft is easily obtained with a high degree of confidence [Arthur D. 
Little, 1966]. Using 0.02 F/ft for the class interval, the frequency 
distribution of gradients was plotted for various XBT traces by signifi- 
cant levels. As expected, the greatest range in gradient values usually 
occurred in the significant level just below the mixed layer where the 
range of temperature values AT was the greatest. 

This method appears to give a good indication of the variability 
of the thermal structure in the depth intervals investigated. Several 
examples of the different types of gradient distributions, using histo- 
grams and a cumulative frequency diagram, are shown in Figures 25 through 
31. The actual analyzed XBT is also shown on the left of each figure. 
Looking at data off San Diego, the first three gradient distribution 
diagrams (Figures 25, 26 and 27) show fairly similar cumulative fre- 
quency distributions. The next two (Figures 28 and 29) show nearly 
identical distributions for water inside the Gulf Stream. The next two 
gradient distributions (Figures 30 and 31) show water in the Gulf Stream's 
frontal region. These last two distributions are not identical for many 
reasons, but mainly because the positive gradient region in the second 


example is below the 100 to 200 foot interval shown in the gradient 
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distribution plot. These results indicate that the small scale structure 
may be water mass related. However, this could not be fully examined in 


the time alloted for this thesis but is recommended for further 


investigation. 
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V. APPLICATION OF NUMERICAL DERIVATIVE MODEL 


This chapter describes the application of the numerical derivative 


model or method (NDM) to several cases of real and assumed data. 


A. BASIC CASES 

Graphs are used to show the features of the NDM. For each station, 
temperature versus depth (T-Z) and temperature gradient versus depth 
(G-Z) plots were drawn by a CalComp plotter using the computer program 
with the DRAW subroutine in Appendix B. The scales are constant for all 
plots. The grid squares were drawn for 100-ft depth increments from the 
surface to a depth of 900 ft. On the T-Z plots (reporduced digitized XBT 
data), temperature was scaled at l1OF per grid line. On the G-Z plots, the 
Eradient was scaled to 0.1F/ft per grid interval, with the axis offset so 
that two-thirds of the plot was allocated to negative gradients. This 
scaling was chosen so that each hundredth of an inch, the resolving 
power of the CalComp plotter, represented one foot of depth and 0.OIF 
of temperature or 0.001F/ft of temperature gradient. The plots included 
here have been reduced in size, but the relative proportions remain the 
same. 

1. Linear Gradients 

For a temperature distribution with linear thermal gradients, 

the G-Z plot is characterized by steps indicating a constant gradient 
during a temperature change. For this thesis the strength of a gradient 
is defined as the absolute value of the gradient (in F/ft units) and is 
described as strong when greater than 0.3F/ft or weak when less than 


0.05F/ft. The term signed strength is defined as the true gradient 
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value and will thus distinguish absolute and true gradient values. The 
magnitude of a gradient is defined as the vertical extent of the gradient 
measured in feet. The product of magnitude and signed strength gives the 
total temperature change for the extent of the gradient. 

a. Layer Structure 

Figure 32 shows a two-layer temperature structure model with 
a T-Z plot on the left and a G-Z plot on the right. The strength of the 
upper layer is 0.05F/ft with a magnitude of 200 ft. The lower layer is 
weaker (0.02F/ft compared to 0.05F/ft) than the upper step but the magni- 
Gude 16> ereaeer (G00 ft to; 200 fe) 

For the multi-layered model (Figure 33) several linear 
gradient steps are shown. In the gradient region between 70 and 90 ft# 
the gradient is strong and goes off scale (Figure 33b). 

b. Step Structure 

A step structure is characterized by thin strong gradient 
regions. The model (Figure 34) shows clearly that the four steps bet- 
ween 60 and 100 ft are stronger than any other steps shown. These four 
steps also have the smallest magnitude (five ft) in the figure. The 
steps at 200, 300, and 400 ft all have the same strength but increase 
in magnitude with depth. 

The features of the T-Z trace (Figure 34a) are easily 


described by the G-Z model. The product of signed strength and magnitude 


gives the total temperature change from the surface temperature to a 


given depth for a step-structure temperature proftigier The temperature 
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Figure 32. 
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Figure 33. Multilayer Model. 
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Figure 34. Step Structure Model. 
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2. Continuous Curving Gradients 


Since finite values are used for both depth and temperature, a 
truly continuous curve cannot be drawn. The spacing by the CalComp 
plotter of every hundredth of an inch can cause variations of this 
size. Computed values of depth, temperature and temperature gradient 
were provided as inputs to a CalComp plotter at every hundredth of an 
inch to create a quasi-continuous curve, at least to a resolution of the 
plotter. It is this quasi-continuous field which accounts for the 0.01 
inch steps encountered in curves drawn by the plotter for any truly con- 
tinuous curve. This effect is also present whenever the plotter is 
drawing any line not directly in the X or Y axis directions. Accordingly, 
any steps shown of this size could be considered as possible adjustments 


made by the plotter or data to adjust to the hundredth of an inch grid. 





An example of a Gaussian profile (Figure 35) shows a normal 


distribution of temperature with depth for the function, like 


iS 
4 fico 
Se 


Z 
T(z) = fe e dx (8) 


where: x = a dummy variable 


80 = sea surface temperature. 
Thé G-Z curve highlights the changes in the T-Z curve. The increasemene 
gradient strength is very clear, reaching a maximum 200 feet, whereas 


the inflection point of the T-Z plot is difficult to find by visual 


inspection. The summation of signed strength values at every foot of 
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Figure 35. Gaussian Model. 
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depth from the surface down gives a total temperature change which is 
analogous to the product of magnitude and signed strength with all 
magnitudes set equal to one. With this summation and a given sea surface 
temperature (SST), a continuous temperature profile can be reconstructed. 
— 
Using an exponential expression for the temperature distri- 
bution 


-nZ/1 
m2 jes Comes) 


(9) 
where: 

60 = SST, 
the T-Z (Figure 36a) and G-Z curves (Figure 36b) appear as compliments 


of each other. The derivative of T(Z) with respect to Z gives: 


-11Z/ 1000 
dZ 50 


which is the equation of the G-Z plot. 
3. Superpositions 
In the example of layers superimposed on an exponential model 
(Figure 37), the conditions for the exponential (9) hold for the depth 
intervals; surface to 40 ft, 220 te 340 ft and 460 to 800 ft yeu 
various layers at the remaining depths. Isothermal layers exist below 
the weak surface exponential gradient layer and also in the 140 to 180 ft 
depth layer. From the T-Z plot, the layers at 340 and 400 ft are not as 
evident as they are in the G-Z plot. 
4. Inversions 
A model with two inversions and several linear steps (Figure 38) 
shows the positive gradient regions clearly. The second deeper inversion 


is clearly a weak gradient, but it is still discernible in Figure 38b. 


68 





Temperature (fF) Gradient (F/ft) 
50 60 -.3 -0.2 -0.1 0 





(2) = 60 - 20e°™ 71000 | 


| 100 
| 200 


| 300 


Zz 
(ft) 


| a 
| 500 
| : 600 


| 700 
a 

| 

| 800 


Figure 36. Exponential Model. 
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Figure 3/7. 
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Exponential Model with Layers Superimposed. 
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Figure 38. Multi-Step Model With Two Inversions 
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With inversions, summing the signed strength values for each foot gives 
the temperature change from the surface and adding this value to SST 


yields the actual temperature at any depth. 


B. OBSERVATIONAL CASES 
Several different oceanic regions were selected to establish fea- 
tures of this model. The oceanic regions explored were =alon= gem. Atlantic 


Coast,,in the Andaman Sea, Eastern North Pacific, and the Gulf Stream. o 


) 
1. Atlantic Coastal Region 


Stations in Massachusetts Bay occupied during December 1969 


were used. The stations were within a radius of 13 miles. Water depths 


were between 280 and 420 feet. “Lack of strong gradients in all seven 
‘traces shows that the water was probably well mixed (Figures 39 through 


42). A positively signed strength in the upper 300 feet was common in 





the traces. Inversions (Figure 40) are quite clear in both the T-Z and 
G-Z plots. A portion of the actual XBT trace is shown for three of the 
stations. In these three cases, there is excellent visible correlation 
between the actual XBT trace and the reconstituted T-Z plots which shows 


that MICROXBT can resolve very small temperature features. In no case 





was there an isothermal layer greater than 140 feet thick. 
2. Andaman Sea 
In the Andaman Sea, near the Straits of Maalacca, a series of 
15 XBT stations from the 18 to 20 July 1969 period were used. The 
stations were within a radius of 12 miles. Water depths varied between 
290 and 380 feet. This oceanic region could be considered somewhat 


Similar to the Atlantic Coastal region above. 
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gure 39. T-Z and G-Z Plots for Atlantic Coastal Region 6 and 7 December 1969. 
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Figure 40. T-Z and G-2 Plots for Atlantic Coastal Region 7 December 1969. 
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ure 41. T-Z and G-Z Plots for Atlantic Coastal Region 9 and 10 December 1969. 
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T-Z and G-Z Plots for Atlantic Coastal Region 11 December 1969. 
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Figure 42, 
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The T-Z plots (Figures 43 through 50) contain a mixed layer to 
200 or more feet except for one observation at 1200Z on the 19th (Figure 
48b). The mixed layer transients are marked with spikes in the G-Z plots. 
Figure 50 shows transients with the greatest magnitude in the upper 100 
feet. The seasonal thermocline is characterized by two steps (Figure 
44, for example). The G-Z plots in most cases show this second step, 
averaging at a depth about 300 feet, with good resolution. Some strong 
gradients are found in the 200 to 300 feet depth region in every G-Z 
plot except in Figure 48b. This figure appears to indicate flow over a 
hump or have erroneous depth values. The fall-rate of the XBT probe, 
appearing erroneously in the T-Z plot, indicates qa possible malfunction — 
of the probe or the recorder for this one cast. 

3. Pacific Oceanic Region 

A series of XBT traces, taken over a week period in October 1968, 
in the deep waters (over 100 nm off the West Coast near San Diego, 
California) were used for this Pecervaevont A composite T-Z plot of 
the actual XBT temperature-depth traces is shown in Figure 51. This 


region is sub-tropical oceanic [Tully, 1964]. Plots of T-Z and G-Z 





(Figures 52 through 65) are presented from the surface to 900 feet. The 
seasonal thermocline is evident in all plots. The strong thermocline is 
clearly shown. Some transients are notable. Few inversion layers were 


present, and, generally, they were found only in the near-surface layer. 


The mixed layer contained numerous temperature fluctuations in 


most cases. Some inversions, in particular the one at 0600Z on 14 October 
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Figure 43. T-Z and G-Z Plots for Andaman Sea Region for 18 July 1969. 
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ure 44. T-Z and G-Z Plots for Andaman Sea Region for 18 July 1969. 
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Figure 45. T-Z and G-Z Plots for Andaman Sea Region for 18 and 19 July 1969, 
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ure 46. T-Z and G-Z Plots for Andaman Sea Region for 19 July 1969. 
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Figure 47. T-Z and G-2 Plots for Andaman Sea Region for 19 July 1969. 
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gure 48. T-Z and G-Z Plots for Andaman Sea Region for 19 July 1969. 
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Figure 49. T-A and G-Z Plots for Andaman Sea Region for 19 July 1969. 
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Figure 50. T-Z and G-Z Plots for Andaman Sea Region for 20 July 1969. 
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Composite T-Z Plot from Actual XBT Traces Down to 900 feet 


for Region 100 nm off 


Figure 51. 
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sure 52. T-Z and G<Z Plots for Eastern North Pacific for 1800Z 10 October 1968. 
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Figure 53. T-Z and G-Z Plots for Eastern North Pacific for 0600Z 11 October 1 
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T-Z and G-Z Plots for Eastern North Pacific for 00002 12 October 1968. 
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Figure 56. T-Z and G-Z Plots for Eastern North Pacific for 1800Z 12 October 1968. 
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Figure 5/7. 


2 








Z(210- Jett 


131058 ~~ Jd600 322 3h 1202?H 


Figure 58. T-Z - 
and G-Z Plots for Eastern North Pacific for 06002 13 October 1968 
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T-2 and G-Z Plots for Eastern North Pacific for 1200Z 13 October 196 
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igure 60. T-Z and G-Z Plots for Eastern North Pacific for 1800Z 13 October 1968. 
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Figure 61. T-Z and G-Z Plots for Eastern North Pacific for 0600Z 14 October 19€ 
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igure 62. T-Z and G-Z Plots for Eastern North Pacific for 12002 14 October 1968. 
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Figure 63. T-2 and G-Z Plots for Eastern North Pacific for 1800Z 14 October 1968 
98 


ea 
iu™ 
co 








| 
L 

I 

- ALM 
ry 


i 

i 

i 
| 
f 
| 
I 





x10 


151088 0000 3240N 120401 


igure 64. T-Z and G-Z Plots for Eastern North Pacific for 0000Z 15 October 1968. 
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Figure 65. T-Z and G-Z Plots for Eastern North Pacific for 0600Z 16 October 1968 
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1968 at 120 feet was clear on both the bathythermogram and the T-Z plot. 
The steps are quite irregular on the T-Z plot but not as irregular as on 
the G-Z plot. 
4. Gulf Stream Region 
Thirteen stations, during the period 11-13 August 1969, were 
analyzed for the Gulf Stream region (Figures 66 through 80). The ship's 
track crossed the Gulf Stream boundary several times. 
a. Interior of Gulf Stream 
For the interior of the Gulf Stream (Figures 70 and 71), 
the mixed layer is distinct and a typical seasonal and main thermocline 
are shown. A "boot" in the 100 foot layer is clearly shown by positive 
gradient values in several of the figures and the actual XBT traces. 
b. Frontal Region 
Strong gradients are the rule for this region and they are 
quite distinct in both the T-Z and G-Z plots. The contrast between 


the frontal and interior regions of the Gulf Stream are quite apparent. 
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T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. 


Figure 66. 
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T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. 
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Figure 67. 
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Figure 68. T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. Actual XBT T-Z 
Trace Inset. 
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Figure 69. 
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Figure 70. T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. Actual XBT T-Z 
Trace Inset. 
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Figure 71. 
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Figure 72. 
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T-Z and “G-Z Profile for Transects across Northern Boundary of the Gulf Stream. 
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T-2 and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. 
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Pigure 74, T-2 and G-Z Profile for Transects across Northern, Boundary of the Gulf Stream. 
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Figure 75. T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. 
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Figure 76. T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. 
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Figure 7/7. T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. 
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Figure 78> ¥-2 and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. 
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Figure 79. T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. 


2(x10° it 


PS 


J SoUsee 
Soe 





Z(x107) ft | 
Figure 80. T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. 
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VI. CONCLUSIONS 


The results give a meaSure of the microthermal structure obtained 


with an XBT by the NDM. 


A. PHYSICAL SIGNIFICANCE 

Microthermal structure statistics have not been adequately explored. 
This thesis has proposed an easy method to observe the fine-scale thermal 
structure. Further, this proposed objective analysis scheme does not use 


the T-Z values but rather the G-Z values. It is the G-Z distribution 


which is important. ‘The T-2 can always be reconstructed if one additional 


1. Internal Waves - Vaisala Frequency 


The Brunt-Vaisala frequency is the maximum frequency of oscil- 
lation of free internal waves. It may be calculated from the density 
gradient, or from the parameters which define density gradient. In 
isohaline water, the NDM gives a direct measure of density (or stability). 
In such waters it could be used to set an upper limit on the maximum 
frequency of oscillation of internal waves. Further, a time series 
study would clearly illustrate the amplitudes likely to be encountered 
since the sheets indicated on the G-Z diagram are readily identified. 

A correlation study of G-Z, Brunt-Vaisala versus depth diagrams would 
be useful. 
2. Oceanic Fronts 

This difficult region to analyze microscopically could possible 

be better explored by using the NDM for the microthermal processes. The 


NDM explores the continuously changing conditions in frontal regions. 
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Frontal regions are usually indicated by extremely complex G-Z plots 
which illustrate the extreme gradients, both positive and negative, 
encountered. Based on the admittedly small number of cases examined, 
there appears to be a tendency of the layers to maintain their character- 
istics and slide by, over or under each other rather than mix. 

3. Sound Propagation 

The effect of microthermal structure on sound propagation 
require better statistical knowledge of the oceans. The NDM gives 
emphasis to the strongest gradient regions. 

The microthermal effects on index or refraction variations 
(density fluctuations) must be determined. Using the NDM with salinity 
variations (from an STD cast) and the NDM with sound velocity variations 
(from an SVTP cast) taken from the two instruments simultaneously, cross- 
correlation techniques could shed some light on the unresolved micro- 
processes in the ocean. In particular, what are the exact physical and 
theoretical correlations among the salinity, temperature, and sound 


velocity. 


B. SUGGESTED STUDIES 
Possible studies for the future have suggested themselves while 
developing this thesis. 
1. Mathematical Models 
Various mathematical models attempted could be explored further 
because each method had some merit. Suggested models are limited by 
ENE IEE ty. 
a. Analysis Methods 
The microthermal structure recorded by XBTs can be looked at 


as a "roughness" factor, the noise level of the water structure. Further 
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work is required to determine such a factor for practical operational use. 
This roughness factor should be computed as the XBT trace is being 
analyzed at FNWC by a program like DIGITBT. This factor could be one 
number to describe objectively the microthermal structure variations on 
the G-Z trace. 

Variations in the NDM should be used with the variables that 
have been arbitrarily set. For example, the depth interval of AT/f£t 
at one foot should be varied to find the effects a larger or smaller 
depth interval would have on the G-Z profiles. As the steps become 
smaller, the exact derivative is approached. Larger steps should smooth 
the spikes shown in the G-Z figures. In some cases, this smoothed G-Z 
profile might be described with curve fitting techniques as suggested by 
Boston [1966] and Grosfils [1968]. 

MICROXBT's digitizing interval of 0.01 inches should be 
changed to 0.005 inches and the program modified to use this interval. 
This would yield more detailed thermal structure than MICROXBT. Obviously, 
one cannot make the interval finer than the noise level of the sensor or 
recording instrument. 

b. Physical Processes 

The G-Z profiles from the NDM emphasize the changing physical 
features which are obscure in the T-Z profiles. By using G-Z profiles 
for parameters other than temperature, clarity in determining distri- 
butions of chemical parameters would be similarly emphasized. Limitations 
in use of the NDM on other chemical properties should be further 
investigated. 

Water mass identification is possible with G-Z profiles. 


Additional work is required in this area. 
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2 Experimental Work 


A microthermal sheet probe could be designed for investigating the 
chemical and physical processes inside the sheet and inside the layers. 
A comparison of actual physical property differences should be made. 

A modified expendable temperature gradient (XTG) meter could be 
designed using two thermistors to obtain the temperature differences in 


the probe. The standard XBT system equipment already installed could be | 





changed to give analog G-Z profiles to the operational user without any 


delay. 
a. Time Series 
The XTG meter could be used to construct G-Z variability at 
one location by conducting a time series analysis. This needs to be 


done in order to obtain a roughness factor, an objective microstructure 
measuring parameter. 

b. Transect 

An XTG meter could be used in transects across various water 
masses to obtain the definite characteristics for the water masses. The 
XTG could be used to follow interesting features on the G-Z plots. More 
important, however, is that the spacial characteristics of microthermal 
structure could easily be determined. 
3. Handling of Data 

When the microthermal structure roughness factor is determined, 
this factor should be made an integral part of the records for that XBT 
trace. This information should be included as part of the ID data. 

Aboard Navy ships the AN/UNQ-7 tape recorder could be used to 
obtain magnetic tape records of the analog XBT or XTG signals. These 


tapes could be submitted with the traces to FNWC or NPS for further analys 
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4. Naval Applications 


Results from microthermal studies show that high frequency sonar 
performance as influenced by the microstructure receive the greatest 
microstructure effects. More knowledge of the microstructure is needed 
before the torpedo and mine-hunting sonar performance could be improved. 
With a measure of the microstructure, the performance of a sonar system 
could be predetermined. 

Tens of thousands of XBT traces are on file at FNWC. New data 
are arriving daily. This thesis is an attempt to show how this data can 
be used both operationally and scientifically. Obviously it represents 
only a beginning in the way of microstructure analysis and operational 


Werlization. 
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APPENDIX A 


MICROXBT FORTRAN IV PROGRAM FOR CDC 6500 
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MICROTHERMAL STRUCTURE ANALYSIS FORTRAN SOURCE PROGRAM 


FOR IBM 360/67 
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JOB CONTROL LANGUAGE (JCL) FOR READING FNWC TAPES 
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//(Job name) 
ey EXEC 
//FORT.SYSIN 


/* 
//G6.FTO6FOO1 
//G6.FTO4FO01 
// 

// 
//GH.FTO4F002 
Vy 

// 

//GO.SYSIN 


JOB (Acct. no., Job no., Section), ‘(Your name)' Note 
FORTCLGP , REGION=100K, TIME . G#=2 Note 
DD * 


FORTRAN Source peck 


Note 
DD SYS@#UT=A ,SPACE(CYL, 6) Note 
DD UNIT=2400-1,V#L=SER=NPSXXX , [ABEL=(,NL), Note 
DCB=(DEN=1 , RECFM=F , LRECL=136 , BLKSIZE=136 ,TRTCH=ET), '' 
DIS P=(@LD, KEEP) Note 
DD AFT=FTO4F001 ,V@#L=SER=NPSXXX , LABEL(2 ,NL), Note 
DCB=(DEN=1 ,RECFM=F , LRECL=136 ,BLKSIZE=136,TRTCH=ET), '"' 
DISP=(@LD, KEEP) Note 
DD * sy Note 


DATA Cards, if any 


/* (Orange Card) 


Notes 


1. Standard JOB, green card, for program submission. 


2. If not plotting or card punching, the FO@RTICLGP can be changed 


to FORTCIG. 


3. Must be a white card. 


4. Optional card for additional storage for output. 


5. The XXX following NPS is the NPS tape number. These three cards 
are for reading the first file on the input tape. These cards 
assign input unit 4 for the source programs READ statements. 
Further note that the record length is 136 characters, all of 
which may be read. This is the standard record length from the 


Cpc 3600. 


6. For a third file to be read, these three cards must be repeated 
with FOO2 changed to FOO3 and IABEL(2,NL) changed to IABEL(3 ,NL) 
and so forth for additional files. Also the UNIT=2400-1 should 
be changed to AFT=FTO4FO01 for all files after the first file. 
If only one file is to be read, then these three cards are not 


required. 


7. Not required if there are no input data cards. 
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//(Job name) 
//DUMPTAPE 
//SYSPRINT 
//SYSUT1 

// 

//SYSUT2 
//SYSIN 


J@B (Acct. no., Job no., Section), ‘(Your name)' Note 
EXEC PGM=IEBPTPCH 

DD SYS#UT=A 

DD UNIT=2400-1,DISB=@LD, LABEL=(,NL) , V@LUME=SER=NPSXX , Note 
DCB=(DEN=1, RECFM=F , LRECL=136 ,BLKS IZE=136 , TRTCH=ET ) Note 
DD SYS#UT=A 

DD *,DCB=BLKSIZE=80 

PRINT MAXFLDS=1 ,ST@PAFT=1000 Note 
RECORD FIELDS=(120) Note 


/* (Orange Card) 


Notes 


1. Standard J@B, green card, for program submission. 


2. If more than one file is on the tape then separate two cards for 
each file identified by file number by the file number in the IABEL 
(X,NL) where S is the file number, which is implied 1 when omitted. 


3. The STOPAFT limits the number of lines of printed output, in this 
example 1000 lines is used before printing is terminated. 


4. The field is recorded at the maximum amount that this program will 


dump which 


is 120 characters. If there are any characters between 


120 and 133 in the records then this routine will not dump those 


enaracters. 
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